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Abstract 

We examine the prospects of using two alternative and complementary ways to explore the 
regions that are favored by global constraints in two simple unified supersymmetric models: 
the CMSSM and the NUHM. First, we consider BR (B s -> which has recently been 

for the first time measured by LHCb. In the CMSSM we show that ultimate, but realistic, 
improvement in the determination of the observable to about 5-10% around the Standard Model 
value would strongly disfavor the A-funncl region, while not affecting much the other favored 
regions. Second, we show that all the favored regions of the CMSSM will be, for the most part, 
sensitive to direct dark matter searches in future one-tonne detectors. A signal at low WIMP 
mass (< 450 GeV) and low spin-independent cross section would then strongly favor the stau 
coannihilation region while a signal at higher WIMP mass (~ 800 GeV to ~ 1.2 TeV) would 
clearly point to the region where the neutralino is higgsino-like with mass ~ 1 TeV. A nearly 
complete experimental testing of the CMSSM over multi-TeV ranges of superpartner masses, 
far beyond the reach of direct SUSY searches at the LHC, can therefore be achievable. In the 
NUHM, in contrast, similar favored regions exist but a sample study reveals that even a precise 
determination of BR(_B S — > fi + fi~) would have a much less constraining power on the model, 
including the A-funnel region. On the other hand, this could allow one to, by detecting in 
one-tonne detectors a signal for 500 GeV < m x < 800 GeV, strongly disfavor the CMSSM. 



*On leave of absence from Department of Physics and Astronomy, University of Sheffield. 
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1 Introduction 



In November 2012 the LHC reached the end of its current data collecting phase with the proton- 
proton beam at yfs = 8 TeV. A huge amount of data was collected, allowing the CMS and AT- 
LAS collaborations to reach an integrated luminosity of around 23/fb each, and LHCb of around 
2.2/fb. The performance of the detectors at the LHC and the effort of the experimental collab- 
orations have been quite spectacular. The past year brought some experimental results whose 
crucial importance cannot be questioned, even though they still require further investigation 
and confirmation with larger amounts of data. 

Most notably, on July 4, 2012, both the CMS and ATLAS collaborations announced a 5a 
discovery of a particle consistent with the Higgs boson predicted by the Standard Model (SM) 
based on the analysis of 4.9/fb of pp collisions at s/s = 7 TeV [HE]. Both collaborations have 
recently updated their results, combining data from the y/s = 7 TeV and y/s — 8 TeV runs. The 
CMS value of the Higgs-like boson mass, 125.8 ± 0.6 GeV [3J, is based on the analysis of the 
data corresponding to integrated luminosities of 5.1/ fb at yfs = 7 TeV and up to 12.2/ fb at 
y/s = 8 TeV in the 77, ZZ, WW, tt and bb decay channels. The ATLAS analysis combined 
approximately 4.8/ fb of data at y/s = 7 TeV with 5.8/ fb of data at y/s = 8 TeV in the same 
five channels, obtaining 125.2 ± 0.7 GeV 0]. 

On November 13, 2012 the LHCb Collaboration reported the first observation of an excess in 
the rare decay B s — > [i + pT [5]. The measured value of the branching ratio, BR (B s — > n + = 
(3. 2j^J) x 10 -9 , is consistent with the value predicted by the SM. This decay has been long 
considered as one of the best probes for new physics, and in particular for low-scale supersym- 
metry (SUSY), since SUSY contributions can be largely enhanced by the sixth power of tan/3, 
the ratio of the vacuum expectation values of the two Higgs doublets (see, e.g., [6j \7\ \8\ l9l ITU l fTT ] 
for some early studies). 

The agreement of the recent measurement with the SM makes it potentially strongly con- 
straining for the allowed parameter space of SUSY models. On the other side, the result still suf- 
fers from substantial experimental uncertainties - its current 2a upper bound is actually weaker 
than the previous 95% confidence level (CL) exclusion limit BR(B S —> < 4.5 x 10~ 9 

obtained earlier by the same collaboration |12) . 

Finally, on the front of direct SUSY searches, the y/s = 7 TeV and y/s — 8 TeV runs have 
significantly improved the limits on the masses of colored superpartners, allowing this way 
both CMS and ATLAS to exclude increasingly larger ranges of parameters of low-energy SUSY 
models. Currently the most constraining 95% CL exclusion limits on the parameter space of 
the Constrained Minimal Supersymmetric Standard Model (CMSSM) comes from the ATLAS 
search for squarks and gluinos with jets and missing transverse energy in the final states, with 
5.8/fb of data at y/s — 8 TeV [T3]. A similar analysis by CMS based on 11.7/fb of data and using 
the kinematical variable ax as a discriminator is slightly less constraining [14] . On the other 
hand, as we will show in this paper, the CMS razor analysis at y/s = 7 TeV with 4.4/fb [TS] can 
be combined with the most recent olt at yfs = 8 TeV [14] to produce a lower bound on the mass 
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parameters of the CMSSM that, in the regions favored by the global constraints, is comparable 
to the current one from ATLAS. 

In our recent global Bayesian analysis of the CMSSM [16] (as well as in several other recent 
global, Bayesian or x 2 -based, analyses [Hl[Il[Ill23ini[12[13[21[lS]),[3 itwas shown or 
reiterated that, when combining through the likelihood function the BR(£> S — > /i + /i~) bound 
from Ref. [T2], the Higgs mass, limits from direct SUSY searches, the relic density of dark matter 
(DM), an excess in the anomalous magnetic moment of the muon (g — 2)^, and other relevant 
constraints, four clearly identifiable regions of the (too, TO1/2) plane, with too and mi/2 denoting 
the scalar and gaugino soft masses, respectively, remain favored by high posterior probability 
for both signs of the Higgs/higgsino mass parameter /z: 

(a) at small Too < 400 GeV and 600 GeV < mi/ 2 < 1000 GeV, where the correct relic 
abundance is obtained via efficient neutralino-stau coannihilation |26| (stau-coannihilation (SC) 
region hereafter). In this region the lightest bino-like neutralino as the lightest SUSY particle 
(LSP) is fairly light, m x < 450 GeV, and so is the lightest stop, ~ 1 TeV, hence the correct 
Higgs mass is achieved due to maximal stop mixing, A|/m| ~ 6; 

(b) at 1 TeV < m\i% < 2 TeV, where the cross-section for neutralino annihilation is enhanced 
by the s-channel resonance of the pseudoscalar A Higgs boson [27] (A-funnel (AF) region), with 
bino-like LSP in the mass range 350 GeV < m x < 700 GeV; 

(c) for too > 3 TeV, too > TO1/2, in a strip of the (too, TO1/2) plane along the border of the 
non-electroweak symmetry-breaking region, where the neutralino remains bino dominated but 
contains a non-negligible higgsino component [2S| 29; (Focus Point/Hyperbolic Branch (FP/HB) 
region) . In the FP /HB region we found |16j a significantly lower posterior probability primarily 
because it was difficult to obtain the correct mass of the Higgs boson. Also, this region is 
in considerable tension with current 90% CL upper bound from XENON100 [30] on the spin- 
independent cross section cr^ 1 of dark matter (DM) scattering off xenon nuclei. 

(d) in the multi-TeV regime (to > 4 TeV, mi/2 > 2 TeV) there is a large region where the 
neutralino LSP is almost purely higgsino- like 20 . Its mass is almost constant, m x « fj, ~ ITeV 
(1TH region hereafter) so that the relic density constraint is easily satisfied, since for such a 
heavy higgsino LSP coannihilation is no longer effective. 

A similar pattern holds also in the Non-Universal Higgs Model (NUHM), although at some- 
what different locations. Specifically, the 1TH region can be found already at much lower mass 
scales, too < 4 TeV and mi/2 < 2 TeV [31 . (For an updated analysis including LHC data, 
see [25] .) 

Clearly, given such large mass scales most of the favored regions will remain beyond the reach 
of direct searches at the LHC. Only part of the SC and a small fraction of the AF and FP/HB 

1 A note of caution is in order regarding a quantitative comparison of different analyses. First, Bayesian posterior 
high probability credible regions and \ 2 confidence regions need not agree as they are based on two different concepts 
of probability. Secondly, even within the same statistical framework, numerical results often strongly depend on the 
values of input parameters used. For instance, the mass of the lightest Higgs boson in SUSY very sensitively depends 
on the exact value of the top quark pole mass, which is different in, e.g., [16] and |23| . 
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regions will be explored. It is therefore interesting to investigate the power of less direct ways 
of experimentally testing those regions, including projected sensitivities, on the most popular 
constrained SUSY models like the CMSSM or the NUHM. In this paper we will investigate two 
such observational venues: future measurements of BR (B s — » ^ + /-0 at the LHC and expected 
reach of direct search for DM through one-tonne detectors. 

Regarding BR (B s — > in [16 it was also shown that the impact of the experimental 

upper bound (at that time) on BK(B S — > was the strongest on the AF region where 

the SUSY contribution to the branching ratio is comparable to the SM one, while the other 
explored regions were less affected. A similar conclusion was reached in [32] . where the impact 
of the new positive measurement of BR(£> S — > was for the first time investigated in the 

framework of the Constrained Next-to-Minimal Supersymmetric SM (CNMSSM), which also 
features similar favored regions when all the constraints are simultaneously taken into account. 
This points to an interesting relation between BR (B s — > /i + fi~) and the relic density constraint 
in the AF region. 

As stated above, because the current experimental uncertainties are relatively large, the 
positive LHCb measurement of BR(B S — > is actually somewhat less constraining for 

models of new physics predicting an enhancement of the observable than the previous exclusion 
bound. On the other hand, the systematic and statistical uncertainties will be greatly reduced 
when a larger amount of data comes, and are expected to ultimately achieve the level of 5%. It is 
therefore interesting to investigate what impact such projected sensitivities of BR (B s — > /i + /i~) 
will have on the favored regions of the CMSSM and the NUHM. 

Our goal is twofold. First, we will present a Bayesian analysis of the current status of 
the CMSSM for a much broader range of input parameters than in [16]. We will apply the 
most recent experimental determinations of relevant input observables, most notably the Higgs 
boson mass and the top quark pole mass, in addition to the recent positive measurement of 
BR (B s — > We will show that, in the context of the CMSSM, the expected substantial 

reduction of the experimental and theoretical uncertainties in BR(i3 s —¥ /i + /i~) will have the 
potential to strongly disfavor basically the whole AF region. Secondly, we will show that the 
expected reach of direct search one-tonne DM detectors will be able to discriminate between the 
remaining two favored regions, the SC and the 1TH regions. On the other hand, we will show 
that, unfortunately, a similar conclusion cannot be reached in the NUHM because of the freedom 
in adjusting the pseudoscalar Higgs mass tua and fi. Still, in both models one should be able 
to distinguish between the SC and the 1TH regions. Furthermore, any DM signal indicative of 
the AF region would strongly disfavor the CMSSM. 

Recently, Ref. [33] analyzed the impact of the present measurement of BR (B s /i + /i~) 
and its future status on random scans of the CMSSM and on the general MSSM (see also [33] 
for another recent analysis of this constraint in the MSSM), showing that a large fraction of 
the points generated would be excluded once the projected uncertainties in BR (B s — > /i + ^ ) 
are considered. Our study is partly overlapping but differs in some important aspects: 1. Our 
analysis of the CMSSM is performed as a global Bayesian scan, where the constraints are applied 
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simultaneously through the likelihood approach (with the exclusion of XENON100, as explained 
later). 2. Our main goal is to focus on the future ability to use BR(P S — > /i + /i~) to disfavor 
high probability regions of models with parameters unified at the scale of grand unification 
(GUT). As a consequence, we do not investigate the general MSSM, alongside to the CMSSM, 
but rather the NUHM model. 3. Unlike in [33], we will also discuss in the detail the implications 
of future direct searches of DM. 

The paper is organized as follows. In Sec. [2] we will demonstrate semi-analytically how 
BR (B s — > /i + /i~) shows a unique discriminating power over the AF region of the CMSSM. 
In Sec. [3J we will describe our scanning methodology, and highlight the implementation of our 
statistical combination of CMS bounds on SUSY masses. In Sec.|4]we will present our numerical 
results and discussion. Finally, we will give our Summary and Conclusions in Sec. [5j 

2 BR (B s fi+fi -) in the MSSM 

In this section we first quickly review the analytic form of BR(P S — » /i + /i~) in the MSSM and 
next analyze its implications for the AF region of the CMSSM and the NUHM. 

The measurement of the branching ratio is a very good probe of new physics, since in the 
SM the decay rate is helicity suppressed, but can get significant contributions in SUSY. 

A general expression for the branching ratio is [351 13S1 ISZ1 13H] 

Gia 2 Mr th n I 4m?, f / 4m 2 \ „ „1 

(i) 

where Mb, and tb, are the B s mass and lifetime, and Pa, Fp and F$ are the axial- vector, 
pseudo-scalar and scalar form factor respectively. In the SM, F$ and Fp are highly suppressed 
by helicity conservation, and the only remaining term in the curly bracket of in Eq. ([lj is 
1-F/iJ 2 ) where Fa can be expressed in terms of the Wilson coefficient Cio, the muon mass m„, 
and the B s decay constant Jb s , Fa — —im^fB^CiQ. The SM value is still a matter of some 
discussion in the literature, as it is marred by large theoretical uncertainty, due mostly to 
the determination of fg s by the lattice QCD groups. Rcf. 39J estimates the CP-averaged 
branching ratio as BR(P S — > A i+ M _ ) SM = (3.23 ± 0.27) x 1CU 9 , while Ref. [30] gives a different 
value, BR (B s -> m + ^) S m = (3.53 ±0.38) x 10" 9 . 

Notice that the theoretical calculation should be rescalcd by the effects of B s — B s oscilla- 
tions [41) in order to be compared with the experimentally measured value. In this study we will 
adopt the value given in J3j5] for the CP-averaged SM branching ratio and, following [41] , take 
BR (B s — > M + M~)sm = 3-5 x 10 -9 for the value rescalcd by the effects of B s — B s oscillations 
(time averaged). We differ here from [33] where (3.87 ± 0.46) x 10 -9 for the latter was used. 

SUSY contributions to BR (B s — > ^ + A*~) become comparable to the SM when Fs and Fp are 
roughly of the same order as Fa- At the leading order (LO), in the framework of minimal flavor 
violation, the dominant SUSY terms in the Wilson coefficients are given by chargino-squark 
terms only and are proportional to tan 3 f3 [37] . 
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Following the calculation and notation given in |37j , one can write for Fp and F$ 



2 

Fs : p — — ~ M Bb jb 3 Cs,p 



where 



Cs,p — T 



tan 3 /3 
— >lo ■ 



(2) 
(3) 



The dominant contributions to J-lo is given by the charginos and squarks in the loop, 

„2 \ / ™2 \ / ^2 

1 .. . Ill' 

Flo — 

m t sin 6y sin 6* t cos # t 




+ (sgn /j,) m x ± cos 6u { V 2M W sin V 



m t cos 0y sin 9 t cos # t 



where we assumed A22 = lc&1^/(Vt&V^) ~ — A33 = —1, and neglected a term in An 

v ub v:j(v tb v* s ) ~ -10- 2 . 

The -Da(a;) are loop functions, 



D 3 (x) 



:ln: 



(5) 



l-x ' 

and cos^y an d sinflj^y are elements of the chargino mixing matrices defined such that 
UM X ±V T = diag(m x ± , m x ±) (See Appendix B of [37] for notation). 

For the purpose of this analysis we shall assume that the experimental uncertainty in 
BR (B s — > will eventually, with about 50/fb of data at yfs — 14TeV, be reduced to about 

5% [42]. We will also assume that the theoretical uncertainty will reach the precision of 5% [43] . 
Hence BR (B s — » fi + fj,~) ■ = (3-50 ± 0.25) x 10~ 9 with both theoretical and experimental un- 
certainties added in quadrature. Further, we will primarily assume that the current SM value 
will be confirmed by experimental measurements from LHC with the above precision, although 
we will discuss some possible deviations. In particular, we shall briefly discuss the case that the 
current LHCb central value is instead confirmed, BR (B s — > M + M~) pr oj = (3-20 ± 0.23) x 10~ 9 
and the case when the assumed ultimate error will be twice as large. 



2.1 Application to the CMSSM 

In the CMSSM, and more generally in unified SUSY models, Eqs. ([l])-([5]) can be greatly sim- 
plified thanks to relations between the different sparticles. 

The first and the third line in the right-hand side of Eq. Q are always opposite in sign and 
cancel each other out to a good approximation. Since in the SC and AF regions the neutralino 
is strongly bino-dominated, it follows that m ± ~ M2 and m ± ~ ^1, so that sinOu smOy — 
and the fourth line in Eq. Q is dominant. In the FP/HB and the 1TH regions, where the 
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Figure 1: The dependence of BR(5 S -)■ on tan/3 in the CMSSM. Solid blue line: \i > 0; 

dashed red line: \i < 0. Solid horizontal lines: current la error on BR(i3 s — > dashed 



horizontal lines: projected error. In boldface the values of tan /3 typical of each region, (a) SC 



region, (b) FP/HB region, (c) 1TH region, (d) AF region 



roles of m ± and m ± are interchanged, cos Ojj cos 9y — so that the second term is thus 
dominant. In all of the favored regions, by remembering that sm9 t cos9t ~ m t A t / (m~^ — nv-^) 
and that, for moderate to large tan/3, a change in sgn/i implies a change in the sign of sm8jj, 
as tan#;y oc l/(cos/3M 2 + sin/3/i), one can recast Eq. Q as 

m?A t 

■Flo - -^l>z— 2", (6) 

m~ — ra- 
ti t 2 

where T>^ is given by differences of D3 functions and is in general of order 0.1-0.3. 

In our numerical analysis we will use full calculations to higher order, given by the latest 
numerical codes, but we can use the above approximation to show how a projected better 
determination of BR (B s — > fJ, + /i~ ) can affect the four regions favored by the correct relic density, 
as explained above. 

For choices of parameters typical of the SC region, J-"lo is the largest, driven up by large 
values of //, ri > 1000 GeV, and by maximal A t /MsusY, which gives T> 3 ~ 0.3. On the other 
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hand, as can be seen from Eq. (|3]), the branching ratio is suppressed by moderate tan/3 values, 
tan/3 ~ 5 — 30 typical for the SC region. The tan/3 dependence of BR (B s —> for a point 

representative of the SC region (mo = 226 GeV, = 827 GeV, Aq = —1375 GeV) is shown in 
Fig. l jaj The solid blue line gives the case /j, > and the dashed red line the case [i < 0. The 



lines are shown in boldface for tan/3 values characteristic of the SC region. The solid horizontal 
lines give the current la theoretical and experimental uncertainty on the measurement added in 
quadrature. The dashed horizontal lines denote our estimated la projected uncertainties added 
in quadrature. 

In the FP/HB region and, at large mass parameters, in the 1TH region, .Flo is the smallest, 
since A t /MgugY is minimal and m ± f=a fj, < 1 TeV. Moreover, the branching ratio is suppressed 
by large ttia, even if tan/3 can assume a wide range of values. We show in Fig. |l||^b)| the 
tan/3 dependence of BR(B S — > /x + ^t~) for a point representative of the FP/HB region (mo = 
3447 GeV, m,i/ 2 = 866 GeV, A Q — 730 GeV), and in Fig. ] 'c) the same for a point representative 



of the ITeV higgsino region (m = 7989 GeV, m 1/2 = 2854 GeV, A Q = -767 GeV). The color 



code is the same as in Fig. ] 'a) 



Finally, but most importantly for the purpose of this paper, in the AF region tan/3 has to be 
large in order to yield the correct Q x h 2 , as we will explain in the following subsection. Thus, the 
measured value of BR {B s — > fi + /i~) becomes important in constraining the parameter space (/x 
and ttia are comparable to the SC region, while stop mixing is not as large, so that T>3 ~ 0.15- 
0.2). In Fig. |l|[d)| we again indicate with thick lines the ranges of tan/3, for both signs of /i, 
which give Q x h 2 within la (theoretical + experimental uncertainties added in quadrature) of 
the central value. The difference in the allowed values of tan /3 for different sgn /i is a feature of 
the AF region, and was already observed in [16] . We will explain this in the next subsection. 
Notice also that, for [i < the calculated value of BK(B S —> fi + fi~) is more than la below the 
SM value, since the form factors Fp and Fa undergo destructive interference, and one is left 
with a small value of F$. 

2.2 The A-funnel region vs BR(B S 

The AF region is particularly sensitive to the determination of BR (B s — > /i + /i~) because at the 
LO the relic density there depends mainly on the same parameters, tua and tan/3; see Eq. 

AF annihilation occurs when the mass of the pseudoscalar A is close 2m x , so the lightest 
neutralino can efficiently annihilate into SM fcrmions (mostly 6-quarks) through the s-channcl 
exchange of A. For an almost pure bino LSP (higgsino components of order 10~ 2 or less) 
and tan /3 in the range 20-60 (for lower values the channel XX Zh becomes dominant) one 
obtains [57] 

n x h 2 * 3 x 10 ; 2W3 / S , (7) 



where 



const tan 2 /3 

m 2 (4 - m^/m 2 ) 2 + (r^m^/m 2 ) 5 



(8) 
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Figure 2: Green bands: the regions of the CMSSM where vtla = 2m x for — lOTeV < Aq < lOTeV. 
Light gray bands: the regions where av for bino-like DM gives the correct Q x h 2 through A- resonant 



annihilation. Dark gray: the neutralino is not the LSP. (a) fi > 0. (b) fi < 



with the A width being, for [i > (fi < 0), Ta ~ 1-3 (2.0) x 1CT 5 m^tan 2 ^. The constant 
in Eq. ^ depends moderately on kinematical factors, on the neutralino composition, and on 
tan/3. For masses given in GeV, its value is ~ 10~ 25 cm 3 GeV 2 /s. The correct relic density is 
generally achieved for a difference \rriA — 2m x | not exceeding 100 GeV [33] . 

In the CMSSM with bino-likc DM, vha can in principle be close to 2m x for wide ranges of Too 
and TO-i/2, if the value of tan f3 is properly adjusted: the mass of a bino-like LSP is approximately 
given by m x sa 0.44 TOx/2 while vha = 0.72 m 1 / 2 X f ^ ta 5 n ° p ) > where f(x) is a monotonically 
increasing function of x, and f(l) — 1. As a consequence, for large tan/3, with increasing m 1 / 2 
2m x increases faster than uia- Hence, in order to get the resonance for larger m!/ 2 , one needs 
to assume smaller tan /3. However, this does not mean that the correct relic density can always 
be obtained since, as Eq. ^ shows, even for tua — 2m x the cross section becomes suppressed 
with increasing neutralino mass for any given tan j3. 

In Fig. [2] the green bands show the regions of the (mo, mi/ 2 ) plane over which the condition 
mA = 2m x is satisfied for fixed tan/3 and — lOTeV < Aq < lOTeV. Figure 



shows the case 



[i > and Fig. 1 b) ^ < 0. Note that, this is achieved when, for the same tan/3, m 1 / 2 is slightly 
smaller for negative y, than for positive fi. The reason lies in a one-loop tadpole contribution 
to the effective potential of the model [45]. The corrections due to sfermions, charginos and 
neutralinos explicitly depend on sgn /i, leading to positive (negative) contribution to m 2 A for 
H > (/j, < 0). Therefore, when the other parameters of the model are left unchanged, itia 
is slightly smaller for negative /i. For the same reason, values of tan (3 larger than 50 cannot 
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be obtained for negative ix since they would lead to m\ < and no electro-weak symmetry 
breaking. 

Using Eqs. (|7j) and Q, one can now calculate the ranges of that for a given tan/3 

would allow to obtain the correct relic density (within la of the experimental central value) 



when nriA = 2m x . We show them in Figs, 1 '&) and £ b) as gray horizontal stripes. One can 



see that the AF region of the CMSSM is confined to a relatively small part of the (mo, 
plane, where the green and the grey bands of the same tan/? intersect. This is also true in the 
more realistic case where and 2m x are within 100 GeV from one another, as is confirmed by 
numerical scans. As a consequence, tan 8 is also constrained in the AF region: it can take values 
in the range 48-55 for positive fj,, and 38-50 for negative fi. This is the ranges we highlighted 
m boldface in Fig. fljfd)] 

One can also see in Fig. |l|[d)| that the values of mj and tan 8 typical of the AF region are 
the ones that show most tension with the current measurement of BR (B s — > /i + /i~). Moreover, 
it is clear that future, more precise, measurements of the branching ratio will have the potential 
to exclude a broad range of the (mo, TO1/2) parameter space corresponding to the AF region. 

In the NUHM the situation is quite different. The additional soft mass parameters in the 
Higgs sector m_y u and mjy d can be traded, through conditions of electroweak symmetry breaking 
(see, e.g., [31]), for and /i: 

„ m 2 „ — m 2 M tan 2 8 1 „ 

A 4 = ~ -Ml, (9) 



tan 2 8 - 1 2 
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m A = m H d + m H u + , (10) 



which can be adjusted to satisfy the resonance condition, independently of tan 8, for much wider 
ranges of both m and m 1 / 2 - Therefore, in the NUHM the AF region giving the correct relic 
density is not as well localized and occupies a wide part of the parameter space. Also tan 8 is 
now allowed to assume a much wider range of input values, which is crucial from the point of 
view of satisfying the BR(B S —> n + fi~) constraint. We will come back to this point later. 



3 Scanning Methodology and Constraints 

In order to examine the impact of the most recent constraints, including BR(_B S — > /i + /i _ ), on 
the parameter space of the CMSSM and the NUHM we use the Bayesian approach. Wc follow 
the procedure outlined in detail in Refs. [13 EH [TB]. Our goal is to map out the 68% and 
95% credible regions of p(m\d), the posterior probability density function (pdf), given by Bayes' 
theorem, 

= p(d\€(m))n(m) 
PW 

where p(d\£(m)) = C is the likelihood function, which describes the probability of obtaining 
the data d given the computed value of some observable £( m )> which is a function of the 
model's parameters m. C incorporates the information about the constraints, as well as their 
experimental and theoretical uncertainties. Prior probability 7r(m) encodes assumed range and 
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Measurement 


Mean or Range 


Error: (Exp., Th.) 


Distribution 


Ref. 


Combination of: 










CMS razor 4.4/fb , y/s = 7TeV 


See text 


See text 


Poisson 




CMS a T 11.7/fb , ^ = 8TeV 


See text 


See text 


Poisson 


M 


m h by CMS 


125.8 GeV 


0.6 GeV, 3 GeV 


Gaussian 





n x h 2 


0.1120 


0.0056, 10% 


Gaussian 


EH 


S(g~2)l VSY xW 10 


28.7 


8.0, 1.0 


Gaussian 




BR (B -> X s7 ) x 10 4 


3.43 


0.22, 0.21 


Gaussian 




BR (_B„ -> ri/)xl0 4 


1.66 


0.33, 0.38 


Gaussian 


m 


AM Bs 


17.719 ps" 1 


0.043 ps -1 , 2.400PS- 1 


Gaussian 




sin 2 6» oft 


0.23116 


0.00012, 0.00015 


Gaussian 


\m 




80.385 


0.015, 0.015 


Gaussian 


m 


BR (B s -+u+ur) , x 10 9 

\ ' ' / current 


3.2 


+1.5-1.2, 10% (0.32) 


Gaussian 





BR(B S -j-zi+tt") . x 10 9 

V r r / pro.) 


3.5 (3.2*) 


0.18 (0.16*), 5% [0.18 (0.16*)] 


Gaussian 






* We will also consider the case of projected uncertainties around the current measured central value. 



Table 1: The experimental constraints that we apply to constrain model parameters. 



distribution of m. Finally, p(d) is the evidence and is a normalization constant as long as only 
one model is considered, but serves as a comparative measure for different models or scenarios. 

Bayes' theorem provides an efficient and natural procedure for drawing inferences on a subset 
of r specific model parameters (including nuisance parameters), or observables, or a combination 
of both, which we collectively denote by ipi . They can be obtained through marginalization of 
the full posterior pdf, carried out as 

p(A=i,..,r\d) = J p(m\d)d n ' r m, (12) 

where n is the total number of input parameters. An analogous procedure can be performed 
with the observables and with a combination of the model's parameters and observables. 

3.1 Experimental Constraints 

The central object in our analysis is the likelihood function as the place where theoretical 
predictions are compared with experimental data. The constraints that we include in the current 
analysis are listed in Table [TJ As a rule, following the procedure developed earlier [54], we 
implemented positive measurements through a Gaussian likelihood, in which the experimental 
and theoretical uncertainties were added in quadrature. For the Higgs mass, we used the most 
recent CMS determination of its central value and experimental uncertainty, as it is in perfect 
agreement with the determination obtained by ATLAS at the end of the y/s — 8 TeV run. The 
theoretical uncertainty was estimated to be 3 GeV [TBI E!] ■ 

As stated above, for BR (B s — > we considered two cases: 

1. The current measurement at LHCb, for which we adopted a theoretical uncertainty of 
10% of the measured value (see next-to-bottom row in Table [TJ, in agreement with [40] once the 
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uncertainty due to the top pole mass (~ 1%) is subtracted. We do so because in our scans the 
top mass is one of the nuisance parameters and the effect of varying it is included parametrically. 

2. The projected 'best-case' scenario for the determination of BR(_B S — > [i + fi~), where the 
experimental and theoretical uncertainties are both reduced to 5% of the measured value (see 
bottom row in Table [TJ, as explained in Sec.[2j In addition, as a sensitivity test, we considered 
both the case where the measurement will be narrowed down to the time-averaged SM value, 
3.5 x 10 -9 , and the case where the current central LHCb experimental value, 3.2 x 10 -9 , will be 
confirmed by future sensitivities. This second case can in principle improve the fit for the AF 
region in the fi < case, since the branching ratio there assumes values more than la below the 



SM determination (see Fig. 1 [d) and |16|). Finally, we will double the assumed error around 
the SM value, again as a sensitivity test. 

Following the procedure already adopted in our previous papers, we did not include the 
XENON100 upper bound explicitly in the likelihood function. The theory uncertainties are 
very large (up to a factor of 10) and strongly affect the impact of the experimental limit on the 
parameter space. The main source of error (the so-called E^at term [56]) arises from different, 
and in fact partly incompatible, results following from different calculations based on different 
assumptions and methodologies. Such uncertainties do not follow a particular statistical dis- 
tribution, and are not well suited for inclusion in a likelihood function. Moreover, we showed 
in a previous publication [47] that, when smearing out the XENON100 limit with a theoretical 
uncertainty of order ten times the given value of o~p the effect on the posterior is negligible for 
regions of parameter that appear up to one order of magnitude above (and below) the experi- 
mental limit. However, even if we do not include the XENON100 bound in the likelihood, below 
we shall comment on its possible effects on the posterior pdf. 

The likelihood for limits from direct SUSY searches deserves a more detailed explanation, 
which we give in the following subsection. 



3.1.1 Combination of CMS SUSY search limits 

In previous work [35] [TB] we presented a methodology for deriving approximate but accurate 
likelihood functions for two of the direct SUSY searches with all-hadronic final states at CMS: 
otT (Ref. [46 ) and razor (Ref. 16J). Our approximation correctly reproduced the 95% CL 
exclusion bounds of those searches in the (mo, fn\i<i) plane. In [32] we then showed that the 
same procedure for the razor search could be extended to the CNMSSM. 

The likelihood maps were developed through a step-by-step procedure which included gen- 
eration of the SUSY signal at the scattering level with PYTHIA6.4 [57] and a simulation of the 
CMS detector response with PGS4 [SB] to calculate the efficiency once the kinematic cuts were 
applied. The obtained signal yields were finally statistically compared to the publicly available 
observed and background yields of the searches to construct the likelihood map. 

As mentioned in the Introduction, the most constraining limit for the CMSSM presently 
comes from the ATLAS search for squarks and gluinos with jets and missing transverse energy 
in the final states, with 5.8/fb of data at y/s = 8TeV |T3J. The recent limits produced by the 
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CMS Collaboration with comparable or larger luminosity P3] are slightly weaker. On the other 
hand, riding on our accurate method for constructing the likelihood function for all-hadronic 
SUSY search limits with the information provided by the CMS Collaboration, we are in a 
position of deriving an approximate statistical combination of the CMS searches at yfs = 7 TeV 
and y/s = 8 TeV. 

We prefer to follow this procedure rather than taking the ATLAS limit as a hard cut (as 
recently done, e.g., in [59j E2]) for one important reason. Most recent analyses of the CMSSM 
have pointed out that the region of parameter space which provides the best fit to the constraints 
(particularly the Higgs mass) is the SC region. Since this is the region directly adjacent to the 
exclusion bounds, accurate modelling of the likelihood function becomes important. 

In what follows we briefly summarize the methodology adopted for the razor in our previous 
papers, since it will be used again here. We then proceed to statistically combining it with the 
most recent CMS olt search to update our exclusion bound. 

Razor 4.4/fb, yfs = 7 TeV 

The CMS razor search, based on 4.4/fb of -Js = 7 TeV data, found no excess of events over 
the SM prediction. In deriving the likelihood map for the razor analysis we followed the CMS 
procedure described in [15] . All accepted events were divided into 38 separate bins in the two- 
dimensional space of the razor variables R 2 and Mr, and the likelihood of observing a certain 
number events in a given bin was defined as a Poisson distribution convolved with a Gaussian 
or log-normal function that would take care of the predicted error on the background yields. 
The details of our analysis can be found in [TB] . 

a T 11.7/fb, yfl = 8 TeV 

The CMS «t search, performed with 11.7/fb of data based on <Js = 8 TeV pp collisions, shows 
no significant deviation from the SM prediction [T5] . In deriving the likelihood map we followed 
closely the CMS procedure and our methodology presented in [IB] H?]. The accepted events 
were divided into 8 separate boxes, according to the number of jets originating from 6-quarks, 
rib = 0, 1, 2, 3 or > 4 and to the number of reconstructed jets per event, 2 < rij < 3 and 
rij > 4. In every box, the events were classified based on the value of the variable Ht, defined 
as the sum of all jets' transverse energies. The likelihood for observing Oi events in the i-th 
bin, given the known number of the expected events Sj, and the number of the expected SM 
background events is given by a Poisson distribution convolved with a Gaussian, to account 
for the predicted error on the background yield. The ranges of Ht in every bin, together with 
the corresponding numbers of the observed events, expected background events, and errors on 
the expected background yield provided by the CMS Collaboration, are given in [60] . 

Since the CMS Collaboration has not provided for this search the official 95% CL exclusion 
bound in the CMSSM, we validated our likelihood map procedure for the ay analysis with the 



official ar 4.98/fb, y/s = 7 TeV contour given in [ST]. In Fig. i '&) we show the comparison 



between our simulation and the official plot. As one can see, we obtained very good agreement. 
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Figure 3: |(a)jThe 68.3% CL (red solid thin), 95.0% CL (red solid thick), and 99.7% CL (red dashed 
thin) exclusion bounds for the CMSSM from our approximation of the ax likelihood (y/s = 7 TeV, 
~ 5/ fb) compared to the orig inal CMS 95% CL exclusion bound (dashed black), [(b)] 95% CL 
exclusion bound (solid red) for the CMSSM from our approximation of the ay likelihood (y/s = 
8 TeV, ~ 12/fb) compared with the bounds obtained for the NUHM, when mu u < m H d (dashed 
blue) and mn d < mn u (dot-dashed black), (c) 95% CL exclusion bound (solid red) for the CMSSM 



from our combination of CMS searches (solid red) compared the current ATLAS bound (dotted 
gray). 
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In Fig. c b) we show our 95% CL contour for the olt 11.7/fb, \J s = 8TeV as a solid red line. 
It is a big advantage of the likelihood map methodology that it allows one to derive likelihood 
functions for SUSY searches even where the official limits are not available. 

Furthermore, we also show in Fig. |3^b)| that the derived exclusion limit can be applied not 
only to the CMSSM, but also to the NUHM. The exclusion bounds obtained for two different 
choices of the parameters mu u and mjj d (shown in dashed blue and dot-dashed black) do not 
differ from the CMSSM one. The reason is that the soft masses of the Higgs sector enter the 
one-loop renormalization group equations of the first two generation quarks only by the terms 
multiplied by the Yukawa couplings, and therefore are strongly suppressed, while the term pro- 
portional to the difference {m 2 H — rrtg ) is multiplied by the factor <? 2 /10 and is also negligible, 
unless the mass difference is very large. The NUHM exclusion limits shown at Fig. |3j[b)| corre- 
spond precisely to the choice of parameters that would maximize the difference \rrifj — m 2 ^^ 
and at the same time remain in agreement with the physicality condition. 

In [T5J we showed that the 95% CL limit based on the 4.4/fb razor search is not affected by 
the change of the sign of parameter ^. The same is true for the NUHM. 

Limit combination procedure and results 

In our approximate combination of the recent SUSY searches by CMS, we used all bins consid- 
ered in the CMS ax 11.7/fb, y/s = 8TeV analysis [13], as well as the ones from the razor 4.4/fb, 
■s/s = 7 TeV [16] . Following the statistical approach of Modified Frequentist Confidence Lev- 
els [52] we assumed that the two searches are statistically independent (since they are based on 
different data sets) and we treated every bin as a statistically independent counting experiment. 
Then the combined likelihood is a product of the likelihoods for the two separate searches. The 
results of such an approximation are presented in Fig. c |cj which shows a comparison of the 



95% CL lines for the 5.8/fb ATLAS search at s/s = 8 TeV (dotted gray) and our combination 
of CMS results described above (solid red) . 



3.2 Scanning tools and parameter ranges 

In this analysis we used the package BayesFITS which calls several external, publicly available 
tools: for sampling it uses MultiNest [53] with evidence tolerance factor set to 0.5, sampling 
efficiency equal to 0.8, and number of live points equal to 4000 (CMSSM) or 10000 (NUHM). 

Mass spectra were computed with SOFTSUSY v3.3.6 [53] and passed via SUSY Les Houches 
Accord format to Superlso v3.3 65J to calculate BR (B — > X s "f) , BR (B s —> /i + /i~), BR (B u —> rv), 
and £( 9 -2)* u& \ AM Bs , sin 2 8 e s and M w are calculated with FeynHiggs [55J . DM observ- 
ables, such as the relic density and direct detection cross sections, are calculated with Mi- 
crOMEGAs 2.4.5 [57]. 

The prior ranges and metric adopted for scanning the CMSSM and nuisance parameters are 
given in Table [2] We only scanned in log priors for the mass parameters, as it was proven in 
many previous studies [55] that flat priors in the CMSSM unduly favor the large-scale regions of 
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CMSSM parameter 


Description 


Prior Range 


Prior Distribution 


m 


Universal scalar mass 


0.1, 20 


Log 


mi/2 


Universal gaugino mass 


0.1, 10 


Log- 


A 


Universal trilinear coupling 


-20, 20 


Linear 


tan/3 


Ratio of Higgs vevs 


3, 62 


Linear 


sgn/i 


Sign of Higgs parameter 


+ 1 or -1 


Fixed 


Nuisance 


Description 


Central value ± std. dev. 


Prior Distribution 


Af t 


Top quark pole mass 


173.5 ± 1.0 


Gaussian 


m b (m 6 )slf 


Bottom quark mass 


4.18 ±0.03 


Gaussian 


a s (M z ) W 


Strong coupling 


0.1184 ±0.0007 


Gaussian 


l/ Qcm (M z ) W 


Reciprocal of electromagnetic coupling 


127.916 ±0.015 


Gaussian 



Table 2: Priors for the parameters of the CMSSM and for the SM nuisance parameters used in our 
scans. Soft masses and Aq are in TeV. Top quark pole mass and bottom quark mass are in GeV. 



the parameter space (volume effect). Moreover, the correlation between BR(i? s — » and 
the AF region, which we expose in this study, becomes unobservable once large values of tua 
become favored by the scan. 

Note that, compared to [15], we significantly extended the ranges of Too, m 1 / 2 and Aq. Wc 
performed our scans for /i > and /i < separately. For negative /i we did not include the 
(g — 2) M constraint, since its only effect would be to worsen the overall fit (see jl6j for a detailed 
discussion of this issue), while the observable is very poorly fit anyway. 

The prior ranges and metric for the NUHM parameters are given in Table [3] We performed 
several scans with different choices of ranges and priors. As we will explain in more detail in the 
next section, we selected the ranges that allowed us to most strongly see the possible correlation 
between BR(B S — > /i + /i~) and the AF region. 

4 Results 

In this section we will present our numerical results. We will first examine the impact of the 
current and the projected determination of BR(_B S — > fi + fi") on the different high probability 
regions of the parameter space of the CMSSM and will discuss ensuing implications for testing 
them. In particular, we will show that the AF region is likely to be basically fully excluded if 
the SM (or else current) value of BR (B s — > /i + /i~) is confirmed with high precision. Next we 
will demonstrate that future one-tonne detectors of dark matter scattering off nuclei will provide 
a crucial complementary way of cross-examining those regions and of potentially exploring the 
favored regions of the CMSSM over very wide ranges of parameters not accessible to direct LHC 
searches for new particles. Next we will apply a similar approach to the NUHM and show that 
the above conclusions in general will not hold. On the other hand, some positive measurements 
of DM signal will have the potential to basically rule out the CMSSM. 
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NUHM parameter 


Description 


Prior Range 


Prior Distribution 


m 


Universal scalar mass 


0.1, 4 (0.1, 20*) 


Log (Linear) 


mi/s 


Universal gaugino mass 


0.1, 4 (0.1, 10) 


Log (Linear) 


A 


Universal trilinear coupling 


-7, 7 (-20, 20) 


Linear 


tan j3 


Ratio of Higgs vevs 


15, 35 (3, 62) 


Linear 


sgn A* 


Sign of Higgs parameter 


+ 1 or -1 


Fixed 


m » u 


GUT-scale soft mass of H u 


0.1, 4 (0.1, 20) 


Linear 


m Hd 


GUT-scale soft mass of Hd 


0.1, 4 (0.1, 20) 


Linear 


Nuisance parameters like in the CMSSM 



* In parentheses we show the ranges for the scans giving the ITH region, see Sec. |4.2| 



Table 3: Priors for the parameters of the NUHM and for the SM nuisance parameters used in our 
scans. Soft masses and Aq are in TeV. Top quark pole mass and bottom quark mass are in GeV. 



4.1 The CMSSM 

In Figs. 4^a) and |^[b)| we plot 68% and 95% credibility regions of a two-dimensional (2D) 
marginalized posterior pdf (henceforth called posterior for brevity) in the (mo, m.1/2) plane and 
in the (Aq, tan/3) plane, respectively, for /i > 0. In Figs, ^[c) and \4j[d)\ we show the same for 
/1 < 0, but without 5 (g — 2) , as mentioned earlier. The figures give an update and a significant 
extension of the results presented in our previous CMSSM analysis [16] by an inclusion of the 
new positive measurement of BR (B s — > /i + /i~) (instead of an upper limit) and by significantly 
extending mass parameter ranges; compare Table [2] In the ranges overlapping with those in 
the previous study (0.1 TeV < m < 4 TeV, 0.1 TeV < m 1/2 < 2 TeV, -7 TeV < A < 7 TeV 
and 3 < tan/3 < 62) the figures basically reproduce the same features, with the main three 
regions of high posterior favored primarily by the DM relic density and the Higgs mass, and 
also by the other constraints. First, the SC region shows up at small too just above the LHC 
(CMS, and similarly for ATLAS) exclusion linej^] The posterior features a 68% credibility and 
the best-fit point is located there thanks to a very good fit to the Higgs mass, and a value of 
BR (B s — > in agreement with the experiment (at not too large tan/3). Next, the AF region 

can be seen at the 2a credibility level for 1 TeV < m < 4 TeV and 1.2 TeV < m 1/2 < 2 TeV, 
although a much smaller la 'island' at smaller mo is also present. Finally, the FP/HB region 
appears only as a 95% credibility island at mo « 4 TeV due to the fact that it is more difficult 
there to produce the correct Higgs mass. (See [H] for a detailed discussion, and also [32] where 
we discussed in detail the CMSSM limit of the CNMSSM, and adopted the same updated values 
of experimental constraints as in this study.) 

As a side remark, we note that in [16] the best-fit point was located in the AF region]^] With 
the new improved fit the best-fit point is now found in the SC region - this is due to the updated 



2 The SC strip can be narrowed down by applying limits on long-lived charged particles to staus |25| but this will 

not significantly change the results presented here. 

3 It was also emphasized there that the location of the best-fit point in the CMSSM is very sensitive to exact values 

of input parameters, approximations used, etc. 
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(a) the (mo, m^) plane of the CMSSM for \i > 0, 



(b) ) the (A), tan/3) plane for \x > 0, (c) the (mo, m^ 2 ) plane for fi < 0, and (d) ) the (vlo, tan/3) 
plane for \x < 0, constrained by the experiments listed in Table |TJ with the exclusion of 5 (g — 2)^ 
for \i < 0. The 68% credible regions are shown in dark blue, and the 95% credible regions in light 
blue. The dashed red line shows the CMS combined 95% CL exclusion bound. 



(somewhat increased) value of the top pole mass which made it easier to obtain a 126 GeV Higgs 
mass in the SC region, also in the CNMSSM, as we discussed in detail in [32] . 

In the case of /i < (but without <5 (g — 2) ) the AF region is much less prominent than 
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for /i > 0, although still visible in Fig. ^c) at 95% credibility. Likewise the FP/HB region has 
shrunk considerably, while the SC remained fairly stable. 

Going to larger mo and m^, beyond those considered in |16j . the main new feature in 
Figs, ^a) and ^[c) is the appearance of a large 68% posterior region ranging from around 5 
to 12TcV in mo [M] where the LSP is an almost purely higgsino-like neutralino with mass 
m x s» fj, ~ 1 TeV (the 1TH region). The correct Higgs mass is also easily obtained there due 
to large -Msusy while all other constraints, including BR(£> S — > (except at large tan/3), 

reproduce basically the SM value there. In fact, for /i < with S (g — 2) dropped from the list 
of constraints, the best-fit point has now moved up to the 1TH region since no other constraint 
favors lower MgusY- 

Notice that Figs. |^[b)| and |^[d)| show that the parameters Aq and tan /3 are now less con- 
strained than in |16j . This is a consequence of extending the scanned ranges of tuq and TO1/2 
to much larger values. The large higgsino DM region corresponds to two large 68% credible 
regions, where tan/3 assumes values in the range 30-55, and Aq can take very large negative 
and positive values. The tree level value of rrijj is often positive, in which case electro-weak 
symmetry breaking (at the tree level) is not achieved. To overcome this, large and negative 
one-loop contributions to m 2 H , which are proportional to |^4t|tan/3, are needed. Therefore, at 
smaller tan/3 larger \A t \ (and hence \Aq\) are favored - a tendency that becomes weaker as tan/3 
grows. Hence, the posterior features a 'gap' for small Aq, which narrows down with increasing 
tan/3 (0 < Aq < 5 TeV implies \A t \ <C Aq for most choices of the other parameters). Values of 
tan/3 < 25, on the other hand, are not favored in this region given the prior ranges considered 
in this study. 



In Figs. I 'a) and I b) we show the marginalized 2D posterior pdf in the (mo, TO1/2) 
plane for fi > and fi < 0, respectively, for the scans where we adopted projected future 
theoretical and experimental uncertainties (added in quadrature) for BR (B s — > fi + fi~), i.e., 
BR(B S —> M + M~)proj = (3-5 ± 0.25) x 10~ 9 . The AF region does not appear in the high pos- 
terior anymore, even at 95% credibility, while the other regions basically do not change. This 
was to be expected in light of the discussion presented in Sec. [2] The argument is valid for 
both signs of /1: when fi > (fi < 0) BR(£? S — > /i + /i~) in the AF region assumes much larger 
(smaller) values than the ones favored by the projected uncertainties, as shown qualitatively in 
Fig. |]|[d)| The location of the best-fit point is different for /1 < since, like in Fig. ^'c) the 
constraint from <5 (g — 2) has not been included in the likelihood function. 

The results shown in Fig. [5] are quite insensitive to the projected uncertainties assumed for 
BR (B s — > ^i + /i _ ). The shape of the posterior pdf does not change even if they are doubled, 
which we have checked numerically. We also point out that the situation will not be different 
if in the future the uncertainties are narrowed around the currently measured central value 
(3.2 x 10 -9 ) instead of the SM value. We also checked this numerically, finding no significant 
difference in the resulting posterior. 

To highlight the fact that in the CMSSM a more precise determination of BR (B s — > /i + /i~) 
can lead to an almost complete exclusion of the AF region, we show in Fig. [6] scatter plots 
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Figure 5: Marginalized 2D posterior pdf in the (mo, ffiift) plane of the CMSSM constrained by 
the experiments listed in Table[l]with projected uncertainties for BR(1? S — > (a) [i > 0, |(b) 



[i < 0. The 68% credible regions are shown in dark blue, and the 95% credible regions in light blue. 
The dashed red line shows the CMS combined 95% CL exclusion bound. 
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Figure 6: Scatter plot of the points in the (mo, m^j) plane of the CMSSM for 
fi < 0, satisfying Q x h 2 at 2a (pink squares), Q x h 2 + BR(i? s — > /x + /i~) at 2a (blue circles), and 
Q, x h 2 at 2a and \rriA — 2m x | < 100 GeV (green triangles). 
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of the points of our chains restricted to the low-mass regions. In Fig. (. '&) fj, > 0, whereas 
in Fig. ( [b) /i < 0. Pink squares mark the points for which the relic density constraint is 
satisfied at 2a (theoretical + experimental errors added in quadrature); blue circles represent 
the subset of these points for which BR(B S — » /i + /i~) is satisfied at projected 2a (with total 
a = 0.25 x 10~ 9 ); green triangles mark the subset of these points that belong to the AF region 
(\rriA — 2m x | < 100 GeV). One can see a good spacial separation between the blue and green 
points, which is a reflection of the tension of the AF region with the BR (B s — » constraint. 

The mass scales typical for the AF region are so high that most of it will remain beyond direct 
reach of the LHC. (For example, with 300/fb at 14TeV, CMS will probe m 1/2 < 1.3 - 1.4 TeV 
at m < ITeV [69].) Likewise, the FP/HB region will also be only partially probed at the 
LHC, while the 1TH region will remain completely beyond direct collider reach. As we have 
demonstrated, the projected precision in the determination of BR(_B S — > /i + /i~) will have the 
power to potentially rule out the AF region, but not the other ones. 

Fortunately, an expected ultimate sensitivity of DM searches in deep underground detectors 
will provide a crucial complementary, and partly overlapping, venue of testing all the high 
posterior probability regions of the CMSSM. In Figs. 1 a) and |7|[b)| we show a 2D posterior 



we can 



in the (m x , a^ 1 ) plane for /i > and /i < 0, respectively. Starting from Fig. 1 [&) 
clearly identify the four high posterior probability regions, each with a characteristic LSP mass 
range and er^ 1 . The SC region (appearing only at 95% credibility but featuring the best-fit 
point) corresponds to fairly low m x (< 450 GeV) and typically the lowest or 1 . Next to it, with 
somewhat larger m x and a^ lies the AF region (also at 2a credibility), which in turn is very 
well separated from the big 68% credibility region of ~ 1 TeV higgsino LSP. Finally, the spin- 
independent cross section in the FP/HB region featuring a mixed bino- higgsino neutralino (a 
horizontal branch at erf 1 ~ 3 X 10 8 pb) already shows tension with the current 90% CL upper 
bound from XENON100 [30] ■ However, as demonstrated in |47 and mentioned above, this region 
is probably not yet firmly excluded due to large theoretical and astrophysical uncertainties. The 
remaining three regions are currently below the XENON100 exclusion line but will be almost 
entirely probed by future one-tonne detectors, as a projected sensitivity line for XENON-IT [70] 
indicates. Note that, in the absence in the likelihood function of any constraint to favor the 
SC or AF regions, the broad ranges of the CMSSM input parameter assumed for our scans the 
posterior strongly favor the 1TH region which presents the vast majority of points (the volume 
effect) even with the log prior on mo and mj/j, although at 95% of total posterior probability, 
the other regions are also present. We also note that in random scans one can find points with 
reasonably good x 2 (Sx 2 < 12) lying also beyond those favored regions. We illustrate this by 
superimposing on the posterior a distribution of samples uniformly selected from our nested 
sampling chain (blue dots). 

For comparison with the situation at present, Figs, i jcj and |^dj| show the same posterior 
in the case where the future projected uncertainties on BR(£? S — » are assumed; in other 

words BR(_B S — » /i + /i~) is assumed to be basically reproducing the SM value. For [i > the 
AF region is now gone and there remain essentially two testable regions: the 1TH region, which 
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Figure 7: Marginalized 2D posterior pdf in the (m x , cr^ 1 ) plane for the CMSSM constrained by 
the experiments listed in Table [I] |(a)| /U > 0, current uncertainties in BR(1? S — > /i + /^~), [(b)] ji < 0, 
current uncertainties in BR(i? s — > fi + fi~), (c) fi > 0, projected uncertainties in BR(1? S — > n + fi~), 

The dashed red line show the 90% CL 



and (d) fi < 0, projected uncertainties in BR (B s — > [i + 
exclusion bound by XENON100 (not included in the likelihood), and the dashed gray line the 
projected sensitivity for XENON- IT. A distribution of samples uniformly selected from our nested 
sampling chain is superimposed. 
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Figure 8: Marginalized 2D posterior pdf in the (BR(i? s — > dp 1 ) plane for the CMSSM con- 

strained by the experiments listed in Table [TJ |(a)| /^t > 0, current uncertainties in BR (B s — > 
(b) /x < 0, current uncertainties in BR(B S ft + /i~), (c) jj, > 0, projected uncertainties in 
BR (B s — > fi + /j,~), and (d) fi < 0, projected uncertainties in BR (B s — > /i + fi~). The dashed 
red vertical lines show the current |(a)| and |(b)| and projected |(c)| and |(d)| uncertainties on 
BR (B s — > at la. A distribution of samples uniformly selected from our nested sampling 

chain is superimposed. 
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should be basically fully reachable by future DM searches only, and the SC region, testable also 
in part by direct searches at the LHC. Furthermore, they are so widely separated in the plane 
that a detection of a DM signal, even with poor initial determination of both m x and up 51 , would 
have the power to discriminate between them. Furthermore, for /j, < the CMSSM predicts 
that only the higgsino region will be reachable by one-tonne detectors, while in the SC region a 
well known cancellation of two terms reduces up 51 to hopelessly low values. This actually gives 
one a chance, even if somewhat indirect one, to additionally determine the sign of \x since any 
DM measurement indicative of the SC region would most likely favor the positive sign of /i. 

Some of the points made above are recast in a somewhat different way in Fig. [8] where we plot 
2D posterior regions in the (BR(_B S — > fi + fi~) , Up 1 ) plane assuming the current (upper panels) 
and projected (lower panels) determination of BR (B s — > as indicated with vertical bars 

showing the combined (theory + experimental) errors. At present the AF region lies (for both 
signs of fx) mostly (at 95% credibility level) beyond the current ler experimental lines, and 
clearly not yet firmly excluded. However, after the projected uncertainties on BR (B s — > /i + /i~) 
are assumed, as shown in Figs. £j[c) and <^d)| only the two testable regions mentioned above 



survive, the 1TH region corresponding to larger up 51 and the SC region at the borderline of 
XENON-IT reach (fi > 0) or below it (fi < 0). 

4.2 The NUHM 

We have demonstrated above that projected but realistic sensitivities of BR(B S —> /i + /i~) will 
have the discriminating power to basically rule out the AF region in the CMSSM. Furthermore, 
future one-tonne detectors of DM will reach down to values of dp 51 such that either a signal 
in one of the two remaining high posterior probability regions is detected, or the CMSSM will 
basically be ruled out over very wide ranges of its parameters (with the exception of the SC 
region at negative /i), thus reaching far above the direct sparticle mass reach at the LHC. On 
the other hand, by detecting a DM signal at low m x the sign of /i could potentially also be 
determined. 

The question arises whether such rather strong statements extend beyond the CMSSM. 
Unfortunately, it is easy to see that this is not the case already in the NUHM, which is one of 
the simplest extensions of the CMSSM. As mentioned above, in the NUHM, one can choose tua 
and ji as the additional two free parameters; see Eqs. ([9])-( 10 1. These are precisely the quantities 



that played the crucial role in the CMSSM where they were, however, tightly constrained. On 
the other hand, we will show that in the NUHM one predicts some signatures for DM searches 
in one-tonne detectors that are absent in the CMSSM - this could provide the way for ruling 
out the latter model over multi-TeV ranges of mass parameters. 

Since the enlarged parameter space of the NUHM, with much more freedom in the Higgs 
sector, allows a very good fit to almost all observables (except invariably 5 (g — 2)^), it is very 
time consuming to perform a global Bayesian scan as above for the CMSSM. Additionally one 
has to worry about much stronger prior dependence and volume effect [31] . However, since our 
goal in this paper is to examine the impact of future BR(B S — > a* + aO an d direct DM search 
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Figure 9: Scatter plot of the points in the (a) (mo, m\/2) plane and (b) (m^, 2m x ) plane of the 



NUHM for \i > that satisfy fl x h 2 at 2a (pink squares), fl x h 2 + BR (B s — > ) at 2a (blue 
circles), and VL x h 2 at 2a and \m& — 2m x \ < 100 GeV (green triangles). 



sensitivities, a scan over a much more limited range of priors, given in Table [3j is sufficient to 
provide a counter-example to the conclusions drawn above in the CMSSM. Furthermore, we will 
not need to draw Bayesian high posterior regions to make our point. 



In Fig. £ a) we present, for /i > 0, the distribution of points (pink squares) in the (mo, 
plane for which the value of the relic density does not exceed the central value by more than 2a. 
(Since the relic abundance is a strong constraint with a very small uncertainty, the distribution 
of points determines 95% credibility regions of the 2D pdf to very good accuracy but we don't 
show them here.) In green we show the subset of these points for which the correct relic density 
is obtained through neutralino annihilation via the A-resonance. These points constitute the 
AF region of the NUHM. We also show in blue the subset of the pink point that will additionally 
satisfy the constraint on BR(B S —> n + n~) within the projected 2a error. In Fig. i b) we show 
the same sets of points in the (iua, 2m x ) plane, to highlight the features of the AF region. A 
very similar pattern emerges for /i < 0, hence we do not show it here. 

By comparing these figures with Fig. |6j one can see that, in contrast to the CMSSM, in 
the NUHM a precise determination of BR (B s — > /i + /i _ ) will have no real discriminating power 
over the regions of the (mo, mj/ 2 ) plane, as the points do not show spacial separation. In 
other words, the AF region will remain prominently allowed even if a future determination of 
BR (B s — > /i + /i~) will narrow it down to basically the SM value. 



and 



The same, unfortunately, is true when it comes to future direct detection of DM. In Figs. 1C a) 
b) we show the same points in the (m x , a?-) plane for /i > and /i < 0, respectively. 



It 



One can see that the green points featuring the AF region with good relic density and SM-likc 
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Figure 10: Scatter plot of the points in the (m x , a^) plane of the NUHM for (a) 
fi < 0, that satisfy Q x h 2 at 2a (pink squares), Q x h 2 + BR(i? s 
O x /i 2 at 2a and \rriA — 2m x | < 100 GeV (green triangles). 



\i > 0, and 



(b) 



) at 2a (blue circles), and 



BR (B s — > /i + /x~) cover very wide ranges of both m x and <7p , extending from the SC region 
to the 1TH region, which is also present in the NUHM [31]; for a recent update see [23] ■ One 
needs to remember that the limited prior ranges used here for the NUHM (see Table [3]) do not 
fully reproduce the large 1TH region at m x ~ 1 TeV. Nevertheless, we added to the plots the 
points of the 1TH region that were obtained with extra scans featuring linear priors in all mass 
parameters, and much broader ranges (also given in Table j3J). This is allowed, as long as we 
do not draw any statistical conclusion from the combination of these chains. One can see from 
Fig. [10] that it will be much more challenging to discriminate among the three high probability 
regions: a DM signal detected at smaller m x could be indicative of either the SC or the AF 
region while the same at m x close to 1 TeV could instead imply either the AF or the 1TH region. 
Furthermore, for negative fj, a large number of cases with good dark matter relic density and 
SM-like BR(£? S — s- (and correct Higgs mass, etc) accidental cancellations produce cr^ 1 

well below the reach of even one-tonne detectors. 

On the positive side, there is one class of DM signal measurements that could potentially 
allow one to basically rule out the CMSSM over a whole reasonable range of parameters. The 
NUHM prominently features a wide region of roughly 500 GeV < m x < 800 GeV and often 
within the reach of one-tonne detectors which is absent in the CMSSM (except for a handful of 



cases with relatively poor x 2 ; compare Fig. 7 ^c) ). A detection of a signal in future DM searches 



indicative of this mass range would then provide a strong argument against the CMSSM. 
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5 Summary and Conclusions 



In this paper we have examined the implications from the current and the projected but realistic 
sensitivities of both BR (B s -> ^ + /-0 at the LHC and a^ 1 in direct DM searches on the CMSSM 
and the NUHM. Within the CMSSM we performed an updated global Bayesian analysis of the 
CMSSM, with particular focus on the impact of the recent measurement of BR (B s — > at 
LHCb. We further extended the parameter ranges with respect to our previous analysis of the 
model, and we updated the limits from CMS direct SUSY searches through our likelihood map 
procedure, obtained by simulating the SUSY signal and the detector efficiencies. We showed 
that the same lower bounds apply to the NUHM as well. We confirmed that, in the CMSSM, 
in addition to the previously identified high posterior probability regions of the (mo, mi^) 
plane favored by the global constraints, at previously unexplored large CMSSM mass scales a 
prominent 68% credibility region appears, where the LSP is a nearly pure higgsino with mass 
of about 1 TcV. 

We highlighted a correlation between BR(B S — > fi + n~) and the A- funnel region of the 
CMSSM as the above branching ratio and the annihilation cross section in the AF region both 
primarily depend on the same parameters: rriA and tan/3. In this regard, we showed that the 
AF region of the CMSSM is at present slightly disfavored (95% credibility of the posterior pdf) 
by the first BR(£? S — >■ /i + /i~) measurement, although far from excluded, given the large exper- 
imental (and theoretical) uncertainties. However, with expected future, significantly reduced 
uncertainties (experimental of 5% of the measured value; theoretical of 5% of the SM value) 
this observable alone will have the potential to basically rule out the whole AF region, and 
thus a very broad range of the (m , m 1 / 2 ) plane that will for the most part remain beyond the 
reach of direct sparticle searches at the LHC. Next we showed that DM direct detection search 
sensitivities expected for future one-tonne detectors provide a complementary and strong tool 
to test and discriminate between the remaining two high probability regions of the CMSSM: the 
SC region corresponding to the LSP mass of < 450 GcV (and borderline cr^ 1 ) and the ~ ITeV 
higgsino region with a much wider range of er^ 1 . Note also that for fi < only the latter case 
is, for the most part, detectable. This also implies that a DM signal indicative of the SC region 
would strongly favor the positive sign of /i. 

The NUHM presents, unfortunately, a much less clear cut behavior with respect to the 
interplay of the above observables. In particular this is so because the pseudoscalar mass can be 
treated as a free parameter of the model and can be adjusted, along with the other parameters, 
in different ways to yield a good fit to almost all observables. While high probability regions 
analogous to the CMSSM are also present in the NUHM, and no additional ones, they correspond 
to different ranges of the parameter space. As a result, unlike in the CMSSM, one can easily 
identify the AF region with very SM-like BR(_B S — > /i+^~). Furthermore, m x and cr^ 1 in the 
AF region extend to much wider ranges than in the CMSSM. For this reason, in the NUHM it is 
unlikely to be possible to use future determinations of BR (B s — > and cr^ 1 to convincingly 

rule out the A funnel, which will also remain for the most part beyond the reach of LHC 
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direct SUSY searches. On the other hand, a measurement of a DM signal in the mass range 
500 GeV < m x < 800 GeV would be a strong indication against the CMSSM where such cases 
giving a good fit too all data are absent. 
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